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Abstract The Fresnel-Fizeau effect is a special relativistic 
effect that makes the speed of light dependent on the veloc¬ 
ity of a transparent, moving medium. We present a theoreti¬ 
cal formalism for discussing propagation of electromagnetic 
signals through the moving Earth atmosphere with taking 
into account the Fresnel-Fizeau effect. It provides the rig¬ 
orous relativistic derivation of the atmospheric time delay 
equation in the consensus model of geodetic VLBI observa¬ 
tions which was never published before. The paper confirms 
the atmospheric time delay of the consensus VLBI model 
used in lERS Standards, and provides a firm theoretical basis 
for calculation of even more subtle relativistic corrections. 
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1 Introduction 

Very long baseline interferometry (VLBI) and satellite laser 
ranging (SLR) are the primary geodetic techniques which 
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are used to measure continental and intercontinental base¬ 
lines with unprecedented precision (SoyerTet^ 1998 Schu^ i 
[and Behrendl |2012[ |Collilieux et al.[ |2007| | . The estimated 
site coordinates for each technique define an international 
terrestrial reference frame (ITRF) [Petit and Luzum[ |20I0| 
Chapter 4]. Currently, ITRF-VLBI and ITRF-SLR are com¬ 
pared by making use of the standard Helmert (seven-parameter) 
transformation which provides the translational, rotational, 
and scale differences between the frames [Torge and Muller 


|2012| . Continuous inter-comparison between SLR and VLBI 
frames consistently indicates to the persisting scale offset 
between the two realizations of ITRF amounting to 1 ppb, 
i.e. 6 mm at the Earth’s surface [Altamimi et al.||2()I I [|Schuh| 


and Behrend 20121. This problem has raised the question 


whether the relativistic ’consensus’ model used by the inter¬ 
national Earth rotation service (IFRS) for VLBI data analy¬ 
sis I Petit and Luzum 2010| is entirely self-consistent. 

It motivated a working group of the lAU Commission 52 
’’Relativity in fundamental astronomy” to start re-examination 
and improvement of the consensus model in several direc¬ 
tions regarding a more appropriate modeling of the VLBI 
time delay equation in the part concerning the gravitational 
effects I Kopeikin et al. 2011[[Soffel and Han| 20141. How¬ 
ever, relativity also couples with the propagation of light 
through terrestrial atmosphere causing a tiny but fundamen¬ 
tally important contribution to the atmospheric time delay. 
Relativistic equation for the atmospheric time delay was in¬ 
troduced to the consensus model by M. Eubanks on a pure 
phenomenological basis [Eubanks] 1 1991) but it lacked a solid 
theoretical justification and was never properly modeled from 
scratch starting from the fundamental VLBI time delay equa¬ 
tion applied in the barycentric celestial reference system (BCRS). 
The present paper is to fulfill this important task. 

At first glance the idea to consider the atmospheric de¬ 
lay in the global barycentric coordinates may look strange 
for the atmospheric delay is related to the local physics of 
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the Earth. It is similarly ’’local” as, e.g., the tectonic activity 
or ocean loading. Moreover, the atmospheric delay model 
(e.g. the Saastamoinen zenith delay and all mapping func¬ 
tions) are defined in a non-relativistic way. This is usually in¬ 
terpreted as a relativistically meaningful formulation in the 
OCRS QPetit and Luzum| |2010| |. For this reason, it seems 
natural to consider the tropospheric delay as an additional 
effect applied to the topocentric moments of reception at 
VLBI stations. In this way, the only relativistic effect which 
one may expect, would be associated with the relative ve¬ 
locity of one VLBI station with respect to another. Indeed, 
this is exactly what Eubanks 1 1991) actually assumed in his 
consensus model of the relativistic correction to the atmo¬ 
spheric delay that appears as the last term in equation Q 
below. 

The problem with the naive ’’geocentric derivation” of 
the atmospheric delay is that it is completely detached 1) 
from the propagation model of light in vacuum, and 2) from 
the set of space-time transformations adopted in the con¬ 
sensus model [Kopeikin et al. 2011 [ [Soffel and Han 2014 


Petit and Luzum 2010|. The vacuum part of the model is 


2 The Consensus VLBI Model 

The relativistic consensus model for reduction of very long 
baseline (VLBI) observations was established about two decades 


done, first, in the BCRS, and then, converted to the OCRS 
moments of observation. Merely adding the geocentric at¬ 
mospheric time delay to the vacuum part of the VLBI de¬ 
lay is unsatisfactory from theoretical point of view for we 
cannot be sure that all, perhaps minuscule but otherwise im¬ 
portant relativistic contributions have been properly taken 
into account. Indeed, the Earth moves along with its atmo¬ 
sphere with respect to the BCRS, and radio signal propa¬ 
gates through the moving atmospheric medium in the barycen- 
tric coordinates. It introduces the specific Fresnel-Fizeau ef¬ 
fect in the time of propagation of light through the atmo¬ 
sphere. This effect was not considered in the consensus VLBI 
model and it remained unclear what kind of the atmospheric 
time delay will remain after transformation of the Fresnel- 
Fizeau effect to the geocentric (OCRS) coordinates. The present 
paper gives an exhaustive analysis of this important theoret¬ 
ical problem. 

The rest of the paper is organized as follows. Section 
1^ gives a brief introduction to the consensus model. Sec- 
tion[^discusses the speed of light in a moving medium and 
Fresnel-Fizeau effect. Sectionj^derives the atmospheric time 
delay in the consensus model and proves that the Eubanks 
phenomenological model of the atmospheric delay and its 
coupling with relativity is correct. Discussion of our results 
is given in|^ 

We use c to denote the constant fundamental speed in 
Minkowski space-time that is equal to the speed of light in 
vacuum. The locally-measured speed of light in medium is 
denoted as C(. Three-dimensional spatial vectors are denoted 
with boldface letters, like a, b, etc. 


ago I Eubanks 19911 and, currently, forms the basis of the 


IFRS Standards [ Petit and Luzum||2010| . The model is de¬ 
signed primarily for the analysis of geodetic VLBI obser¬ 
vations of extragalactic objects (quasars) conducted from 
the surface of the Earth. In what follows, we shall consider 
specifically this case. 

The basic equation of the consensus model yields the 
time delay, f 2 - ft, between the geocentric coordinate times 
(TCG) of arrivals of the front of a radio wave from a quasar 
to the first and second VLBI stations on the Earth’s surface 
I Petit and Luzum||2010 eq. 11.10], 
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where, here and everywhere else, the dot between two vec¬ 
tors denote a usual Euclidean dot product, fyj - fy, is the time 
delay of the radio wave as if it propagated in vacuum, K is 
the unit vector from the barycenter of the solar system to 
the source of light in the absence of gravitational or aberra¬ 
tional bending, are the atmospheric propagation delays 
measured in the local inertial frame at the sites of the first 
{i- 1) and second {i - 2) VLBI receivers, tu, = dxildt and 
Xi = Xi{t) are the OCRS velocity and radius vector of the 
/-th receiver respectively. 


The vacuum delay is given by [Petit and Luzum 2010 
eq. 11.9] 
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where zITgray is the differential gravitational time delay |Pe 


tit and Luzum 2010 eq. 11.7], L® is the barycentric BCRS 
velocity of the geocenter, U is the gravitational potential at 
the geocenter neglecting the effects of the Earth’s mass, y 
is the PPN parameter equal to 1 in general relativity, b = 
X 2 {ti) - xi{ti) is the OCRS baseline vector at the time of 
arrival f j. 

The vacuum part of the relativistic time delay equation 
Q is well-established theoretically and its validity has been 
independently checked many times [Kopeikin et al. 20111. 
The present paper is concerned with the calculation of the 
atmospheric time delay and its coupling with relativity that 















































Fresnel-Fizeau effect in VLBI time delay 


3 


is the second and third terms in the right side of Q- The 
current expression in ([T]) was proposed by M. Eubanks who 
used the idea of the Lorentz transformation of the tropo¬ 
spheric delay between the two sites as measured in the proper 
time of the first observer. However, there is no exact theoret¬ 
ical justification/derivation of these terms in literature and 
it remained unclear why such particular form of the atmo¬ 
spheric time delay in Q was chosen. It is desirable to derive 
the atmospheric delay and its coupling with relativity from 
the first principles without relying upon plausible but in all 
other aspects phenomenological arguments given in support 
of these terms (see [Eubanks 1991 part like]. 

A standard approach for calculating the atmospheric time 
delay is non-relativistic I Petit and Luzum 2010 Chapter 9]. 
It is based on assumption that the atmosphere is static with 
the speed of light, ct - cjn, where n is the refraction index of 
the atmosphere at the point under consideration. The refrac¬ 
tion index n depends on pressure, temperature and humid¬ 
ity of the atmosphere as well as frequency of electromag¬ 
netic wave [Murray 1983 §7]. The pressure and temper¬ 
ature change with the altitude making the refraction index 
a function of spatial coordinates n - n{x). This approach 
does not take into account that calculation of the propaga¬ 
tion of radio wave from quasar to VLBI receivers is done 
in the BCRS. Earth moves and rotates with respect to the 
BCRS, and so does the atmosphere. Hence, when one calcu¬ 
lates the propagation of radio wave in BCRS, one has to take 
into account the relativistic effects caused by the moving re¬ 
fractive medium on the propagation of the radio wave. This 
gives rise to the Eresnel-Eizeau drag effect and introduces 
the atmosphere-relativity coupling terms in VLBI time de¬ 
lay equation. Additional relativistic corrections arise from 
the transformation of the atmospheric time delay from the 
BCRS to observer’s local frame [ [Petit and Luzum| [2010[ 
Chapter 10]. The goal of the present paper is to account for 
the relativistic effects in the propagation of electromagnetic 
wave through the Earth’s atmosphere and to derive the at¬ 
mospheric VLBI time delay from the solid theoretical prin¬ 
ciples. 


3 Speed of Light in a Moving Medium and 
Fresnel-Fizeau effect. 

In calculation of the coupling between the atmospheric time 
delay and relativistic terms it will be enough to take into 
account only special relativity as the gravity effects (e.g. the 
Shapiro delay, etc.) are small during the time of propagation 
of the signal in atmosphere and they can be neglected in 
the atmosphere-relativity coupling terms Let us consider 

* The reader should be aware that the Shapiro time delay in vacuum 
due to the gravitational field of the Earth and calculated on the part of 
light trajectory lying in the atmosphere, is sufficiently large and cannot 
be neglected [Petit and Luzum[|2010[ 


a small element of volume of the Earth atmosphere which is 
moving along with the Earth with respect to the BCRS, and 
rotates with the same angular velocity as the Earth. 

Let u be the BCRS velocity of motion of the volume 
element. Velocity of light in the comoving local frame of the 
volume element is cc - C(k, where ct - cjn is the speed of 
light in the rest frame of the element of the atmosphere, and 
k is the unit vector in the direction of propagation of the light 
ray. Velocity of light in the BCRS is given by the special 


relativistic law of the addition of velocities | Kopeikin et al. 
[20TT] eq. 2.176] 
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where c/, = cip, cl is the BCRS speed of light and p is the 
BCRS unit vector in the direction of propagation of light in 
the moving medium. We notice that the unit vector p ^ -K 
- the difference in the directions of the two vectors is propor¬ 
tional to the discrepancy between c and C( or, in other words, 
n-1. It comes also from the gravitational light-bending (and 
from the parallactic effects even if they are negligibly small 
in the case of quasars). Equation Q predicts that the BCRS 
speed of light, cl, depends on the instantaneous value of the 
3-velocity of motion of the element of the medium, u. 

In what follows, it is sufficient to retain in 0 only the 
linear with respect to the velocity u, terms. This is because 
the BCRS velocity, u, is approximately equal to the orbital 
velocity of the Earth, V^, which ratio to the fundamental 
speed is V^/c ^ 10“^, meaning that the quadratic with re¬ 
spect to velocity u corrections to the atmospheric time delay 
are too small to be measured. Expanding numerator of (3 
a Taylor series with respect to the ratio, m/c, simplifies ( 3 

Cl + U 
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Calculating the absolute value of cl, and expanding the nu¬ 
merator of Q with respect to ujc, yields the speed of light 
in the moving medium. 


CL-ce + (k- u) 


c 

1-i 


(5) 


where the residual terms are of the order of (1 - c^lc^){u^ /c^), 
and have been discarded. 

The last term in the right side of Q has been developed 
in 1830 by A. Eresnel, and confirmed in 1851 by A. Eizeau 
in his famous experiment with light travelling through mov¬ 
ing water | Stachelj 2005) . Before the advent of special rela¬ 
tivity, the effect was interpreted as a proof of the existence of 
partial “aether dragging”. This explanation was abandoned 
in 1907 after Max von Laue | Laue| 1907) showed that spe¬ 
cial relativity predicts the result of the Eizeau experiment 
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from the velocity addition theorem without any need for an 
aether. 


4 The Atmospheric Time Delay 


The Earth atmosphere has two main parts - ionosphere and 
troposphere. The ionosphere contains charged particles - elec¬ 
trons, which make the ionospheric index of refraction, « 
dependent on frequency of the electromagnetic signal. For 
this reason, the ionospheric time delay also depends on fre¬ 
quency. This property is used to calibrate the content of the 
electrons and to eliminate the ionospheric time delay from 
the VLBI time delay model. In what follows, we assume that 
it was done and the main contribution to the atmospheric 
time delay comes from the troposphere. Troposphere con¬ 
sists of neutral particles of air, and its index of refraction, n, 
can be considered to a large extent as frequency independent 


I Murray 1983 §7.4]. 


Theoretical calculation of the VLBI time delay begins 
from the consideration of propagation of light from the source 
of light (quasar) to observer in the barycentric coordinates 
{T,X), where T is the coordinate time (TCB) and X are the 
spatial BCRS coordinates. The light ray path consists of two 
parts - propagation in vacuum, and propagation through the 
Earth atmosphere. The fundamental VLBI equation in the 
BCRS can be written down in the following form |Kopeikin[ 
|T990l[Sh^-Saless et al.|[l99Tl 


T2-T^^ --K ■{X2- Xx)+AT^,^^+AT, 


(6) 


where K is the same unit vector as in Ti is the barycen¬ 
tric time of arrival of the light signal to the /-th receiver, 
Xi are the BCRS coordinates of the receiver at the time T,, 
zlTgrav is the same differential gravitational time delay as in 
Q, and 

^Tatm = Tatm2 “ Tatmi ’ 

is the BCRS differential delay due to the passage of the radio 
wave through the Earth atmosphere which moves along with 
the Earth and rotates. Relativistic transformations of (|^ to 
TCG time t yields 

^ITatm 


tl t\ — tv2 fvi 
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where the vacuum time delay, ty^ - ty ^, is given in Q. Next 
step is to calculate and to express it in terms of the 

TCG time delay, that is . Concluding this paragraph, we 
notice that the TCG time f is a coordinate time scale defined 
in a region of space-time in the Earth’s world tube neigh¬ 
borhood. It relates to the proper time ta measured along the 
world-line of an observer A by an ordinary differential equa¬ 
tion that can be found in lERS Standards QPetit and Luzum[ 


2010 Chapter 10, Eq. 10.6]. 


The BCRS atmospheric time delay, T^^jj^., calculated along 
the i-th light-ray path, is 



where where the integration is performed using T as a pa¬ 
rameter along the light ray, cl is the speed of light in the 
moving atmosphere taken from Q, Tai is the TCB time of 
entering the light ray to the atmosphere, T, is the TCB time 
of arrival of the light signal to the i-th receiver. 

Equation 0 gives 


■ n - (k- u){n^ - l) , 


( 10 ) 

Cl c ' ' 

where the index of refraction, n, differs from 1 by a very 
small amount, n-1 ^ 3 x 10“''^ at normal atmospheric pres¬ 
sure and room temperature | |Murray| |T98^ §7.4]. Hence, 
with a sufficient approximation, - 1 = 2(n - 1), and the 
integral 0 is reduced to 


Ti 

Tatm, (n-l)^l-^fc■^lj^f7’. 

Tci 


( 11 ) 


where all quantities in the integrand are taken on the path of 
the i-th light ray. 



Fig. 1 Geometry of VLBI observation. A plane radio wave from 
quasar travels in the direction shown by the unit vector k. Vectors a:i 
and a :2 point out positions of VLBI receivers at the time of arrival of 
the plane front of a radio wave from quasar to each radio dish. Vector 
r point out to the position of the radio signal at time t and is defined in 
{0 Vector 0)0 is the angular velocity of Earth’s diurnal rotation. 


BCRS velocity, u, of the element of troposphere can be 
decomposed in a linear sum of two terms 

M = (12) 

where is the BCRS velocity of the geocenter, and tUatm 
is velocity of the atmosphere on the light ray trajectory with 
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respect to the geocentric celestial reference system (OCRS). 
Equation ( [T2| ) neglects all relativistic terms of the order of 
1 /c^ which are too small for practical measurement in the 
atmospheric time delay model (but they might be interesting 
for more deep theoretical study of the relativistic propaga¬ 
tion of light in moving media in some astrophysical situa¬ 
tions like rapidly rotating pulsar’s magnetosphere). The tro¬ 
posphere follows rotation of the Earth around its axis with 
the angular velocity, a;, so that we can write 

Iftatm = X r , (13) 

where r is the OCRS radius vector of quasar’s radio signal 
in atmosphere. It can be presented as an algebraic sum of 
two radius vectors (see Eig.[T]), 

r-Xi-(k, (14) 


where at, is the OCRS radius vector of the i-th VLBI re¬ 
ceiver, and { is the geocentric spatial distance between the 
receiver and the current position of the radio signal. 

We substitute ( [T3| l, ( [T4l l to ( [T^ , and notice that k ■ (tu® x k) ■ 
0, while tu® X a;,- = tu,- is the OCRS velocity of the i-th re¬ 


ceiver. Then, equation (111 takes on the following form. 


T = T 

. atmj . atmj 


-k-iYs + Wj) 
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where 

Ti 

TatmY f ^n-l)dT , 

T„i 


(15) 


(16) 


is the TCB interval of time corresponding to the time de¬ 
lay calculated in the static (non-moving) atmosphere in the 
BCRS coordinates. This time interval is not directly mea¬ 
surable, and should be transformed to the time delay, 
calculated in the local inertial frame of observer located at 
the site of the i-th receiver. It suffices to account for only 
special relativistic terms in this transformation which yields 


T = t 

. atmi .atmj 


1 + -k-Y<s + Wi) 

C 


(17) 


5 Discussion 


We have presented an exact theoretical derivation of the at¬ 
mospheric time delay in the consensus model with an ac¬ 
counting for the coupling with (special) relativistic terms. 
We have shown that the calculation of the atmospheric time 
delay requires a proper discussion of the Eresnel-Eizeau drag 
effect © in the BCRS and relativistic time transformation 
formula ([T^ of time intervals between two events (the time 
of entering atmosphere and the time of observation) taken 
on the light ray propagating in atmosphere. 


We have found that the criticism of the tropospheric time 
delay equations given in the works by B. Shahid-Saless et al. 
I Shahid-Saless et al. 1991[ and by I. Shapiro (1983, unpub¬ 
lished) given in | Eubanks |1991 section like] is unfounded. 
In fact, both Shahid-Saless and Shapiro referred their tropo¬ 
spheric time delay to the BCRS, which is equivalent to our 
formula ([^. Taking into account the Eresnel-Eizeau effect 
and the special relativistic time transformation reduces 0 
exactly to the expression ([^l given in |Eubanks 1991). 


SLR/LLR analysis is done in the inertial (non-rotating) 
OCRS. Rotational motion of atmosphere causes the Eresnel- 
Eizeau effect in one-way light propagation model. However, 
in SLR/LLR experiments light propagates troposphere twice 
- from the ground station to retro-reflector and back. The 
unit vector, k, entering the Eresnel-Eizeau effect (|^ reverses 
direction while velocity, u of the troposphere remains prac¬ 
tically unchanged during the travel time of the light pulse. 
Eor this reason, the excess time delay caused by the Eresnel- 
Eizeau effect cancels out in SLR/LLR technique which cor¬ 
responds to the currently used light propagation model. Thus, 
we conclude that the atmospheric time delay and/or its cou¬ 
pling with relativity in the SLR/LLR/VLBI geodetic tech¬ 
niques cannot be responsible for the observed scaling offset 
lAltamimi et al.[ |2011| between the two versions of ITRE 
based independently on VLBI and SLR techniques. 


where we have taken into account that the time transforma¬ 
tion has been performed between the two events taken on 
the light ray. 

Equation 0 is plugged into ( [T6] l and the unit vector 
k is replaced with the unit vector K - -k directed from 
the barycenter of the solar system to the source of light. It 
yields. 


T = t 

. atmj .atmj 


l + -K-Ym + w,) 
c 


(18) 


This expression is used in (|^ to calulate /fTatm that enters 
numerator of the second term in the right side of (j^. It re¬ 
sults in the Eubanks phenomenological equation ([^ thus, 
giving a rigorous theoretical justification for its appearance 
in the VLBI consensus model. 


Einally, we emphasize that the present paper was con¬ 
cerned with the derivation of the Eresnel-Eizeau effect and 
its appearance in the consensus model in the linear, with re¬ 
spect to velocity of the medium, terms. It is rather straight¬ 
forward to apply our formalism for calculation the quadratic, 
and higher-order terms, if necessary. In particular, the linear 
Eresnel-Eizeau effect associated with the velocity of the 
orbital motion of the Earth does not appear in the final form 
of the atmospheric time delay equation 0- It might be inter¬ 
esting to check if the quadratic Eresnel-Eizeau effect associ¬ 
ated with the velocity T® of the orbital motion of the Earth, 
would appear in the atmospheric time delay equation. Such 
a dependence is unlikely from relativistic point of view as 
it contradicts to the principle of relativity [Kopeikin et al. 
2011) but its mathematical proof is still required. 
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